Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 15 April 2010 



(MN IMfeK style file v2.2) 



A multiscale approach to environment and its influence on 
the colour distribution of galaxies 

D. J. Wilman 1 ' 4 ^. Zibetti 2 ,T. Budavari 3 

1 Max- Planck- Institut fur Extraterrestrische Physik, Giessenbachstrafie, D-85748 Garching, Germany. 
2 Max- Planck- Institut fur Astronomie, Konigstuhl 17, D- 69 '117 Heidelberg, Germany. 

3 Department of Physics and Astronomy, The Johns Hopkins University, 3701 San Martin Drive, Baltimore, MD 21218, USA 
4 email: dwilman@mpe.mpg.de 



15 April 2010 



ABSTRACT 

We present a multiscale approach to measurements of galaxy density, applied to a 
volume-limited sample constructed from the Sloan Digital Sky Survey Data Release 5 
(SDSS DR5). We populate a rich parameter space by obtaining independent measure- 
ments of density on different scales for each galaxy. This parameterization purely relies 
on observations and avoids the implicit assumptions involved, e.g., in the construction 
of group catalogues. As the first application of this method, we study how the bimodal- 
ity in galaxy colour distribution (specifically u — r) depends on multiscale density. The 
u — r galaxy colour distribution is described as the sum of two gaussians (a red and a 
blue one) with five parameters: the fraction of red galaxies ( f re d) and the position and 
width of the red and blue peaks (/i rG d, Mbiue, c re d and Cbiue)- Galaxies mostly react 
to their smallest scale (< 0.5 Mpc) environments: in denser environments red galaxies 
are more common (larger f re d), redder (larger /i rG d) and with a narrower distribution 
(smaller cr re d), while blue galaxies are redder (larger fibine) but with a broader distri- 
bution (larger dbiuc)- There are residual correlations of f re d and //blue with 0.5 — 1 Mpc 
scale density, which imply that total or partial truncation of star formation can relate 
to a galaxy's environment on these scales. Beyond 1 Mpc (0.5 Mpc for /i rc d) there 
are no positive correlations with density. However f ro d anti- correlates with density on 
> 2 Mpc scales at fixed density on smaller scales, and /Lt re d anti-correlates with density 
on > 1 Mpc scales. We examine these trends qualitatively in the context of the halo 
model, utilizing the proper ties of haloes within which the galaxies are embedded, de- 
rived by lYang et ail (|2007t ) and applied to a group catalogue. This yields an excellent 
description of the trends with multiscale density, including the anti-correlations on 
large scales, which map the region of accretion onto massive haloes. Thus we conclude 
that galaxies become red only once they have been accreted onto haloes of a certain 
mass. The mean colour of red galaxies // re d depends positively only on < 0.5 Mpc scale 
density, which can most easily be explained if correlations of /i re d with environment 
are driven by metallicity via the enrichment history of a galaxy within its subhalo, 
during its epoch of star formation. 

Key words: methods: statistical - galaxies: evolution - galaxies: haloes - galaxies: 
statistics - galaxies: fundamental parameters - galaxies - stellar content 



1 INTRODUCTION 
1.1 Motivation 

The evolution of galaxies is heavily intertwined with the 
growth of the dark matter dominated structure in which 
they live, as baryons react to their local gravitational poten- 
tial. It is therefore one of the most challenging goals of ob- 
servational astronomy to explain the continuously evolving 
galaxy population in the context of a hierarchical universe. 



Most fundamentally perhaps, galaxies exhibit a distinct bi- 
modality of properties which correlates strongly with the 
properties of the embedding potential. This is true whether 
the potential is measured on galaxy scales via properties 
such as luminosity, stellar or dynamical mass, or on larger 
scales via measurement of the galaxies' local environment. 
In particular, galaxies living in deeper potential wells are 
more likely to have formed stars at e arlier times, both in 
terms of the average stellar age (e.g. iThomas et all 120051 : 
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ISmith et al.l I2006T I and the last episode of star formation 
(and thus a lower fraction are continuing to form stars, e.i 



Lewis et alJl2002 



Kauffmann et al 



Baldrv et aDl2004l ; iBalogh et aHl2004al lbl; 
20041 ). Such galaxies are apparently un- 



able to restart star formation, since they have no stable ac- 
cess to a suitable reserv oir of gas which can cool to form 
a sta r-forming disk (e.g. iGallagher et al.l 1 19751 ; ICaon et al.l 
2000). Simultaneously the morphology of these galaxies 
must evolve from the typical rotation-dominated spiral or ir- 
regular morphologies of star-forming galaxies to form a pres- 
sure domi nated elliptica l , or a bulge plus smooth disk lentic- 
ular (e.g. lDressleilll980l ;l Dressier et al. 1 1 19971 ; I Wilman et al.l 
l2009h . The precise nature of these transformations remains 
controversial, despite evidence from simulations and obser- 
vations for various mechanisms which can explain some or 
all of these observations. Galaxy mergers are especially at- 
tractive in the context of a hierarchical Universe, and can ex- 
plain both morphological changes and the rapid exhaustion 
of cold gas, although additional physics may be required to 
keep a galaxy from reforming a gas disk (e. g. Springel et al.l 
120051 ; iHopkins et"aill2009l ; I Johansson et al.ll2009h . Stripping 
or evaporation of cold and/or hot gas associated to an 
infalling galaxy will take place within a dense hot intra- 
cluster medium (e. g. iGunn fc Gottll 19721 ; IChung et al1l2009l; 
lLarson et all Il980l) . although it is unclear whether these 
process es can be important within lower density environ- 
ments ilKawata fc Mulchaevl |2008| ; iMcCarthv et all 120081 ; 
Ijeltema et al.ll2008h . Tidal interactions between galaxies can 
drive gravitational instabilities within galaxies, with gas loss 
or triggering of star formation as likely consequences. Nu- 
merous fast interactions (harrassment) can also drive mor- 
phological evolution (e.g. iMoore et allll999l ). 

To statistically describe the evolutionary path of galax- 
ies, it is necessary not only to take a firm grasp of the rel- 
evant physical processes but also to place robust observa- 
tional constraints on precisely how the galaxy population 
reacts to its environment. It has been customary in this 
field to use neighbouring galaxies as test particles in order 
to 'measure' the local environment of each galaxy. One com- 
putes a 'local density' by simply counting neighbours within 
some fixed radius and velocity range or, almost equivalently, 
by computing the distance to the ~N th ne arest nei g hbour , 
where N is typically chos en to be 5 or 10 (fDresslerl fl980h; 
IBalogh et ail (|2004al 'l. see lCooper et all (|2005l ); lKovac et al l 
((2009J) for a discussion of these and more complex methods). 
As with so many traditional measures, this measurement is 
fairly arbitrary in nature^, and it is not easy to see precisely 
how it relates to the underlying density field of galaxies and 
dark matter. 

The two point correlation function or power spectrum 
of galaxies on the other hand, measure the excess probabil- 
ity over random of finding two galaxies with a given sep- 
aration. The comparison of galaxy subsets selected by e.g. 
colour, provides a measure of the relative bias of these galaxy 
types as a function of the scale of pair separation, diagnosing 
the scale at which one c l ass dominates relative to another. 
(e.g.lNorberg et alJliooH ; iBudavari et~aill2003l ; IZehavi et all 
120051 ). Marked correlation functions weight the galaxies in 



1 N is typically selected to obtain sufficient S/N whilst retaining 
a "local" measurement. 



each pair by a given proper ty, normalizing by t he unweighted 
correlation function (e.g. ISkibba et alj 120061 ). This exam- 
ines how galaxies with large or small values of the chosen 
property are clustered. The bias or marked correlation as a 
function of scale are intimitely related to the details of how 
galaxy properties depend on their environments. 

These statistics simply provide the mean behaviour 
of galaxies by tracing their overdensity on different scales 
which, however, are themselves closely correlated with one 
another. A better approach would be to examine how galaxy 
properties correlate with overdensity on one scale for a fixed 
overdensity on another, independently computed scale. Now 
that such a large, homogeneous sample of galaxies is avail- 
able with the Sloan Digital Sky Survey (SDSS), there have 
been a few attempts at isolating such behaviour on diff erent 
scales. iKauffmann et "all (|2004h . iBlanton et~aH (120061 ) and 
iBlanton fc Berlindl |2007l ) conclude that galaxies only react 
to their environments on < 1 Mpc scales, with no signifi- 
cant residual dependence on larger scale overdensity. This 
co ntradicts the a p parent dependence on large scales found 
by IBalogh et alj (|2004al ). which may have resulted from 
low signal to noise m easurement of density on small scales 
|Blanton et alj|2006h . These works, however, contained cor- 
related measurements of density on both scales. 

1.2 A New Method 

In this paper we approach this problem in a different way. 
Firstly, we parameterize the environment using annular, 
non-overlapping measurements of density on different scales 
to ensure that they are (formally) independent. Secondly, 
we consider the u — r colour of galaxies and utilize the re- 
markable fact that its distribution is bimodal and can be 
modelled as the sum of two gaussians. This relates to the 
bimodality of galaxy properties, since galaxies dominated by 
old stellar populations have re d u — r colours and actively 
star-forming galaxies are blue. iBaldrv et alj (|2004T ) initially 
demonstrated the power of this technique, illustrating the 
low chi-squared values produced with such a fit to u— r SDSS 
colours, and demonstrating that the luminosity dependence 
of galaxy co lours was strong and easily parameterize d using 
this model. IBalogh et all (|2004bl ) and IBaldrv et ail (|2006Tl 
have since extended this method, with a powerful demon- 
stration of the strong and fundamental dependence of the 
fraction of red galaxies on environment (measured using the 
distance to the 4 th and 5 th nearest neighbours) and stellar 
mass. The parameterization of the double gaussian fit (the 
fraction of red galaxies f ro d, and the position and width of 
each peak) provides five quantities which can be interpreted 
in terms of changes in the spectral energy distribution and 
the potential physical trigger, and which can be correlated 
with environment measured on different scales. 

1.3 Expectations from the Halo Model 

Finally, to ease the physical interpretation of our re- 
sults, we would like to know how measurements of den- 
sity map onto the dark matter dominated 'cosmic web'. 
To this goal we rely on the models that describe the 
clumpy distribution of dark matter in the Universe in 
terms of haloes, i.e. the potential wells in which galax- 
ies form, infall, virialize and interact with one another 
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and with the intra-c l uster /group medium (ICM/IGM) (e.g. 
White fc Reel 1 19781 ; ICoorav fc Shethl |2002l ; iDiaferio et at] 



200ll ). The halo model (which applies on both galaxy and 



group/cluster scales) is a simplification of reality, but has 
been hugely successful in describing the d ynamics and 
large scale structur e statistics of galaxies (e.g. lRubin|[l983l : 
iBerlind et al1l2003h . In particular, the two point correlation 
function appears to demonstrate an inflection at ~ 2Mpc 
which is attributed to the switch from sub-halo scales (the 
one halo term) to halo - halo s cales (the two halo term) 
(e.g. IZehavi et al] 12004 120051 ; ICooravl l2006h . Models of 
galaxy formation and evolution pay homage to the halo 
model by assuming galaxies to care about their environment 
only in terms of their embedding dark matter halo (e.g. 
iKauffmann et al.lll993l : ICole et aimOOOh . This puts a nat- 
ural scale on the expected dependence of galaxy properties 
on their environments: they should not care about scales be- 
yond their halo. However simulations also show that the for- 
mation (collapse) time of a halo depends upon the overden- 
sity o n larger scales (e.g. ISheth fc Torme nf l2004l ; iGao et all 
120051 ). In this case, one might expect galaxies to have evolved 
further in a region of the Universe which is more overdense 
on large scales. 

1.4 Structure of the paper 

Section [2] introduces the sample used in this paper, selected 
from the SDSS Data release 5 (DR5). In Section [3] we de- 
scribe in detail the method used to compute the density 
of galaxies on different scales, paying close attention to the 
completeness corrections required to ensure a robust mea- 
surement in the presence of incompleteness in SDSS DR5. 

Multiscale density information for the 73662 galaxies of 
our photometrically complete sample are available for public 
use on request to the authors. 

The method of fitting a double gaussian model is de- 
scribed in Section [4] Our results are presented in Section [5] 
Firstly we examine the simple dependence of the galaxy 
colour distribution on luminosity and density on a single 
scale in Section 15.11 Then, in Section 15.21 we present our 
main results in which the importance of different scales are 
tested independently. This forms the purely observational 
part of the paper. 

In Section [6] we utilize a catalogue of embedding haloes 
(groups) which has been assigned to SDSS galaxies within 
the context of a halo mode l and with a few implicit as- 
sumptions (|Yang et all 120071 1. We examine the dependence 
of halo-based properties such as mass and halo-centric radius 
on the density computed on different scales and compare this 
to the way in which the galaxy colour distribution depends 
upon the same measurements of density. What emerges is 
a remarkably consistent picture, in which the halo model 
is extremely successful at explaining the colours of galax- 
ies: smaller scale density can drive some aspects of colour 
evolution, but star forming galaxies do not seem to feel the 
halo environment prior to infall. This exercise should pro- 
vide the first, qualitative step towards a fully quantative, 
model-independent description of the way in which galax- 
ies trace their environment. Models can then be constrained 
independently of the assumptions that are required in the 
construction of group catalogues or modeling of the correla- 
tion function. The rich parameter space of multiscale den- 



sities will help overcoming degeneracies and difficulties in 
interpreting single aperture measurements of local density. 

In this paper we assume the cosmological parameters: 
(Q. M ,A,h) = (0.3,0.7,0.75). 



2 SAMPLE 

We base our analysis on the Fifth Data Release 
|Adelman-McCarthv et al.ll2007l . DR5) of the Sloan Digital 
Sky Survey (|York et al.ll2000l . SDSS). We define an initial 
sample of galaxies with a valid spectroscopic redshift de- 
termination, which are part of the so-called 'main galaxy 
sample' (i.e. all extended sources with petrosian dereddened 
magnitude brighter than 17.77 in r-band and mean surface 
bri ghtness within half-lig ht radius fi r ^ 23.0 mag arcsec" 2 , 
see lStrauss fc et al.l 120021 ). This sample covers 5713 square 
degrees. We further restrict the selection to the redshift 
range 0.015 ^ z: the lower limit is chosen such that local 
deviations from the Hubble flow negligibly affect distance 
estimate^]. 

For this sample we compute the two photometric quan- 
tities relevant to the present work, namely the rest-frame 
de-reddened Petrosian r-band (absolute) magnitude and 
the rest-frame model colour u — r. De-reddening is based 
on the estimate of fore ground Galactic extinction derived 
bv lSchleeel et alj (| 1998T ) . while the transformation to rest- 
frame q uantities is performed usi ng the software K-CORRECT 
(v4.1.3 iBlanton fc Roweisl l2007h . based on the five-band 
, q , r , i, z spectral energy distribution (see iFukugita et al.l 



1 19961 , for details on the SDSS photometric system). By 
choosing the petrosian magnitude as a base and proxy for 
the total galaxy luminosity we provide that the SDSS selec- 
tion criteria simply propagate in our sample selection and 
binning, and at the same time we avoid redshift dependent 
biases in the estimate of the total lumino sity (this derive s 



from the definition of petrosian magnitude, Petrosianl 



1976) 



2002), 



On the other hand, as noted in IStouehton fc et al.l 
the so-called SDSS model magnitudes are best suited to de- 
rive colours, as they are computed with a weighted aperture 
photometry scheme common to all five bands, which ensures 
results that are independent on the wildly different signal- 
to-noise ratio in each band. The choice of the u — r colour 
among all possible combinations of the five SDSS pass-bands 
a ppears optimal to st udy galaxy bimodality (as shown, e.g., 
in lBaldrv et aU^OO-lT l: the physical reason is that u — r pro- 
vides a long wavelength leverage (therefore it is very sensi- 
tive to the shape of the optical SED, depending on star for- 
mation history, metallicity and dust extinction) across the 
4000A-break (which is the strongest spectral feature at op- 
tical wavelengths and responds mainly to changes in stellar 
population age). 

The upper panel of Figure [T] illustrates the Petrosian 
r-band luminosity M r as a function of redshift z for a ran- 
dom subsample of our galaxies. The spectroscopic selection 
threshold at a luminosity corresponding to apparent magni- 
tude m r = 17.77 is a strong function of redshift, and a weak 
function of SED type at given redshift (k-correction) . This 



2 In doing so we also avoid photometric problems arising with 
galaxies of large apparent size. 
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Figure 1. Sample selection in top: dereddencd Petrosian luminos- 
ity M r vs rcdshift and bottom: k-corrccted dereddencd Petrosian 
magnitude m r vs rcdshift for a random selection of galaxies in 
the sample. Galaxies in the volume-limited sample (z ^ 0.08 and 
M r ^ —20.0) arc highlighted in red. The blue solid line corre- 
sponds to the 95 %ile of the k-correction distribution at the 
magnitude limit of m r = 17.77. The dashed green line is equiva- 
lent, only for the 50 th %ile of the k-correction distribution, whilst 
the dashed magenta line illustrates the magnitude limit (no k- 
correction) . 



is better illustrated in the lower panel where we show the 
k-corrected magnitude not corrected for the distance mod- 
ulus, as a function of redshift. To assess the effect of the 
k-correction on the luminosity limit, the percentiles of the 
k-correction distribution are computed as a function of red- 
shift, and are applied at the magnitude limit. The solid blue 
and dashed green lines show the 95 th and 50 th percentiles 
respectively, which clearly become larger with redshift as the 
k-corrected magnitude limit diverges from the non-corrected 
limit of m r = 17.77 (magenta dashed line). 

To study environment, it is important that the lumi- 
nosity limit for each galaxy and its neighbours (which will 
be used to estimate the local environmental density) is the 
same, regardless of distance. Therefore, a volume-limited 
sample with a constant lum inosity limit as a fun ction of 
redshift is required. Similar to lBalogh et al.l |2004bl ). we find 
that a suitable compromise between volume and depth can 
be achieved by limiting to z ^ 0.08 and M r ^ —20.0. Galax- 
ies meeting these criteria are plotted red in Figure [TJ and 
within the region delineated by the solid red lines in the 
luminosity-redshift plane (upper panel). All galaxies with 
k-corrections less than the 95 th percentile of the distribu- 
tion, and Petrosian r-band luminosities M r ^ —20.0 will 
be potential spectroscopic targets out to z = 0.08, provid- 
ing a measure of environment. The most extreme galaxies 
with the upper 5% of k-corrections may still be missed at 
redshifts very close to z = 0.08. However these will by defi- 
nition be extremely few and with perhaps the least reliable 
k-corrections. 



The volume-limited sample contains 92082 galaxies, 
with a median redshift of 0.066. 



3 MEASUREMENT OF DENSITY ON 
DIFFERENT SCALES 

The goal of this paper is to test the correlation between the 
galaxy colour bimodality and the environmental density on 
different scales, which we measure as follows. A cylindrical 
annulus, centred on each galaxy in the volume-limited sam- 
ple, is used to measure the local projected density of neigh- 
bouring galaxies above the limiting luminosity between an 
inner and outer radius and within a fixed rest-frame velocity 
range. By selecting non-overlapping annuli we can extract 
complimentary information about the overdensity of galax- 
ies on different scales, which will only correlate with each 
other because the underlying density field correlates on dif- 
ferent scales, coupled with the effects of projection and red- 
shift distortions. For an annulus with inner radius and 
outer radius r Q , the surface density E r ,, ro is given by: 



where N° jIo is the completeness-corrected number of 
neighbours brighter than our limit, within the rest-frame 
velocity range ±dv, and with a projected distance from 
the primary galaxy r where n ^ r ^ r . For the pur- 
poses of this paper we choose neighbours brighter than the 
volume-limited sample depth of M r — —20.0 and within 
dv = lOOOfcms -1 of the prima r y galax y. This is equivalent 
to that used by iBalogh et al.l (|2004bl ), sampling ± ~ la 
for massive clusters in velocity space, or a redshift space 
cylinder of length ± ~ 14 Mpc. On large scales this traces 
the non-linear regime for high densities, whilst voids remain 
underdense (~ 1% of galaxies have zero neighbours within 
r = 3 Mpc). 

3.1 Selection and Correction for Completeness 

In order to estimate local galaxy densities properly, selec- 
tion effects must be taken into account and completeness 
corrections applied. In the first place, we request the photo- 
metric coverage of the region used to compute local density 
to be highly complete. We compute the photometric com- 
pleteness within a 3 Mpc circle centred on each galaxy in 
the volume-limited sample, at the redshift of the galaxy. 
Using the same tools that also power the spatial computa- 
tions wit hin the National Vir tual Observatory's Footprint 
Servic«[f| (|Budavari et alj|2007h . we perform the formal alge- 
bra of the relevant spherical shapes on the sky. The 3Mpc 
circle is first intersected with the sky coverage of the DR5 
survey that yields the observed area. Furthermore we have 
to censor problematic regions due to nearby bright stars, 
bleeding, satellite trails and missing fields. The Sloan col- 
laboration provides spherical polygons of these masks that 
one has to exclude by subtracting their shapes from the 
observed area. In addition, we also censor out bad seeing 

3 \http://voservices.ne t/footprint 
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patches that are w orse than 1.7" in the r band. Using the 
Spherical Library (|Budavaril |2009| ) , we compute the exact 
shape of every cell, and analytically calculate its area. The 
completeness is then just the ratio of this area to the full 
circle. 

The resultant cumulative distribution of photometric 
completeness shows two regimes with 80% of the galaxies 
having completeness above 98.645%. This is an acceptable 
loss of area which will only have minimal effects on our den- 
sity estimates. The remaining 20% of galaxies have much 
less complete photometry within 3 Mpc. This is where the 
edge effects of the SDSS coverage become important. We 
keep the 80% of galaxies with better than 98.645% photo- 
metric completeness, for which the residual incompleteness 
can be attributed to regions removed due to bright stars, 
satellite / asteroid trails and imaging defects. This leaves us 
with a sample of 73662 galaxies. 

Local density estimates are obviously affected by the 
spectroscopic incompleteness of the survey, which is likely 
to depend upon environment. In fact, due to fiber collision 
effects (i.e. fibers could not be assigned to targets within 55" 
of each other) spectroscopic observations are biased against 
galaxies in denser regions of the sky. For example, ~ 10% of 
galaxies with £0,0.5 ^ 23 Mpc~ 2 are missing spectra for half 
their potential neighbours. It is possible to correct statisti- 
cally for the spectroscopic incompleteness effects by comput- 
ing a correction factor for the number of neighbours, based 
upon the fraction of target^] with a valid redshift measure- 
ment. 

The spectroscopic completeness correction factor is 
computed on a range of different scales: 



±y T{ ,r 
1 * r: .r„ 



(2) 



where N?. >ro is the total number of photometrically 
identified spectroscopic targets within the annulus and 
Nrj, ro is the number with redshift. This factor corrects for 
the incomplete sampling of galaxies within dense regions, 
up- weighting the number of neighbours in equation Q] such 
that: 



N! 



x N" 



(3) 



where N"" ro is the uncorrected number of redshift space 
neighbours within our radial, velocity and luminosity con- 
straints. 

Figure [2] illustrates how C 0.5 an d Co 1 depend upon 
So, 0.5 and £0,1 ■ We take the 95%ile of the distribution to 
examine the largest corrections in each bin of density. As 
expected, the largest corrections exist at high density due to 
fiber collision effects, although at low densities the quantized 
number of neighbours combined with finite binning leads to 
artificial selection of high completeness galaxies in particular 
bins (striping). 

Galaxies with no known neighbours ( N"" ro = 0) will 
automatically also have no neighbours in the completeness- 
corrected density (£ ri ,r = 0), regardless of the complete- 



ness correction ( ro ). In theory this means a galaxy with 
many neighbours can still have £ ri ,r = if all true neigh- 
bours have no redshift information. However, only ~ 4% of 
galaxies with one known neighbour within 3 Mpc ( N03 = I) 
have completeness correction larger than 50% ( Co, 3 > 1.5). 
Therefore completeness corrections are typically not large 
for galaxies in low density regions, and it is unlikely that 
the signal from galaxies with £ ri ,r = is dominated by 
galaxies which truly have significantly larger densities. 



4 FITTING METHOD 

We parameterize the gal axy colour distributi o n as a double 
gaussia n, as i llustrated bv | Baldrv et aL l (|2004h . lBalorfr et al.l 
j2004bl ) and iBaldrv et all (|2006T ). to examine the depen- 
dence of the galaxy colour distribution on galaxy luminosity 
and environment on different scales. The colour distribution 
$(C) for a given sample of Nt galaxies is thus fit by the 
following function: 



*(co 
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27T<7bluc 



(4) 

where the five relevant parameters are: the fraction of red 
galaxies f rc( j and the mean colours and widths of the red 
and blue peaks fi rc d, Mbiuc, cried and <Tbi uo . 

Our fitting procedure for any gi ven sample i s base d 
upon the one outlined in Section 4 of IBaldrv et alj (|200-4h ■ 
Here we summarize it and describe the details specific to our 
analysis. 

The galaxies are binned in colour, with a bin width of 
0.1 magnitudes. The fit is limited to the range 0.75 ^ u — r ^ 
3.0. 

The double gaussian model is then fit to the binned 
colour distribution using a Levenberg-Marquardt least- 
squares fifl In doing this, each colour bin is weighted by 
the inverse variance, computed as the square of the Poisson 
error plus a softening factor of two [W = N+2 )- Bins con- 
taining no galaxies are weighted equivalently to bins with 
one galaxy (W = |). 

Low number statistics can produce noisy colour dis- 
tributions, and so additional constraints need to be placed 
upon the fits, and initial estimates of the parameters given, 
in order to prevent the fitting algorithm from converging 
to false minima. Most importantly, we only fit samples of 
at least 25 galaxies, which extensive testing shows recovers 
sensible results. Initial estimates o f fj, and a are taken from 
the fit parameterized in Table 1 of lBaldrv et al. I l|2004l) . An 
estimate of f rc d = 0.5 is applied to the brightest bin with es- 
timates for each successive bin of luminosity taking the value 
computed for the previous bin. Once each luminosity bin is 
split into bins of density, the initial parameter estimates 
are set to the values computed for the full bin. The results 
are insensitive to the precise values selected. The blue peak 
is constrained such that /Ubiue cannot be less than Q.2mag 
bluer than the initial estimate of (J, Ic d- This ensures that 
the two peaks remain separate entities, and do not switch 



4 In this paper "Target" refers to a galaxy which meets the SDSS 
magnitude limit of m r = 17.77 



5 This was implemented using the mpcurvefit procedure within 
idl, written by Craig B. Markwardt. 
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Figure 2. 95%ile of the completeness corrections Cg 5 and Cg j as a function of £rj,0.5 ^0.5,1 ■ Larger completeness corrections 
are required at high density, where fiber allocation is difficult. Quantized neighbours lead to "stripes" of high C s in particular lower 
density bins. 



places. Finally, we also constrain the widths (cr Ie d, crt>iuc) to 
be within the range 0.05 ^ a ^ 0.8. This avoids pathologi- 
cal solutions in cases where the colour distribution exhibits 
large noise spikes or unusually large number of galaxies in 
the tails of the distribution. 

In all cases good fits are achieved, yielding reduced 
chi-squared values near or below unity. As discussed by 
iBaldrv et all (|2004l ). the softening factor of two privides an 
approximation to the uncertainties involved for low number 
counts with a Poisson distribution. Inevitably this approxi- 
mation, combined with the finite weighting of bins contain- 
ing no counts, leads to frequent cases of reduced chi-squared 
significantly less than one. In the following analysis errors on 
the f re d are computed purely statistically, and are insensi- 
tive to (but consistent with) the fitting errors. Errors on the 
other parameters are estimated from fitting. In these cases 
it was necessary to multiply the error on each parameter by 
\fixied) since the errors on these parameters are overesti- 
mated for Xred < 1 (propogated from overestimated errors 
on the datapoints). Whilst this is clearly an approximate 
correction, we note that the fit errors on these parameters 
are only applied in Figures [4] and [5] and that our main con- 
clusions are based instead on the correlation with density of 
properties computed in independent bins. 



5 RESULTS 

5.1 Dependence on Luminosity and Density 

A large number of galaxies per luminosity bin are required to 
examine the multiscale density dependence of galaxy prop- 
erties within that bin. Therefore we limit our analysis to the 
luminosity range —21.5 ^ M r —20.0. 

The top panels in Figures [3] [4] and [5] illustrate the de- 
pendence of the five parameter family (f rc d, /i rc d, Mblue, 
fred, Obiuc) on luminosity, M r , as resulting from the fitting 
procedure applied to the volume-limited samp le within this 
lumin osity range. We confirm the results of IBaldrv et al.l 
(|2004 ) indicating that all five parameters are strong func- 
tions of luminosity: lower luminosity galaxies display a lower 



fraction of red galaxies, bluer mean colours for both peaks 
and a broader red peak but narrower blue peak than do 
more luminous galaxies within the range considered. Small 
differences from lBaldrv et al] (|2004 ) result from improved k- 
corrections (a later version of the k-correct code was used). 
There is no measurable dependence of colour ( /x rc d or /ibiue) 
for constant luminosity galaxies on redshift, indicating that 
the k-corrections are consistent across our small redshift 
range. We note that, for the luminosity range sampled by 
the volume-limited sample, the fraction of red galaxies f rc d 
is only weakly dependent on luminosity, while both [i and a 
show strong variation for both peaks. This luminosity depen- 
dence of these properties within the volume-limited sample 
is consistent with that obtained for a larger sample with 
variable redshift limit as a function of galaxy luminosity. 

Residual quantities (A(f re d), A(/i re d), A(/z b iuc), 
A(cr re d), A(crbiue)) are defined for each subsample (sub) as 
follows: 

A(f red ) = fred( M r , Sub) - f rcd ( M r ) (5) 

A(^) = ^(M r ,sub)-KM r ) (6) 

A(cr) = cr( M r , sub) -cr( Mr) (7) 

The middle and lower panels of Figures [3] [4] and [5] show 
how these quantities depend upon luminosity for three bins 
of density on < 1 Mpc scales (£o,i, middle panels) and 
1 — 2 Mpc scales (£1,2, lower panels). These density bins 
contain the upper, middle and lower thirds of the density 
distribution on these scales, and the whole volume-limited 
sample is represented in black. 

The density dependence of the colour distribution on 
both small and large scales is extremely clear in these figures: 

• The fraction of red galaxies is a strong function of envi- 
ronment, with a difference in f re d between the most and least 
dense thirds of the distribution ranging from 0.15 to 0.19 and 
from 0.09 to 0.15 for densities computed on < IMpc and 
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Figure 3. Top: f rc( j as a function of luminosity, M r , for the full volume-limited sample in the luminosity range of interest. Middle: 
Resid uals for different bins of So i- Bottom: and in different bins of larger scale density, Si : 2- Errors are computed using the iGehrelsl 
jl986h approximation, combining Poissonian and Binomial statistical errors. For double gaussian fits to a large number of galaxies this 
approximates the fitting error, whilst for smaller numbers the statistical error is larger. 




Figure 4. Top: /i ro d and cr rcd as a function of luminosity, M r , for the full volume-limited sample in the luminosity range of interest. 
Middle: Residuals for different bins of So 1- Bottom: and in different bins of larger scale density, Si 2- Errors are pure fitting errors. 
Key as in Figure [51 
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1 — 2 Mpc scales, respectively. This is almost independent of 
luminosity in the range —21.5 ^ M r ^ —20.0. 

• The characteristics of the red peak are also highly signif- 
icant functions of environment, although the dynamic range 
is smaller than for the luminosity dependence. The differ- 
ence in the mean colour /i re< j for the most and least dense 
thirds of the distribution on < 1 Mpc scales is in the range 
0.02 to 0.06 mag. In the most overdense regions, in going 
to fainter luminosity A(/x re< j) increases (becomes redder), 
whilst the relative width of the peak A(o- rc d) simultaneously 
decreases. A similar range but with less obvious luminosity 
dependence is observed using densities on 1 — 2 Mpc scales. 

• The blue peak properties Hblue and (Tbiuc are more 
strongly correlated with environment than the equivalent 
red peak properties. A difference in fx^lne of 0.15 to 0.18 
mag between the most and least dense bins on < 1 Mpc 
scales is present at all luminosities. The difference is still 
significant but slightly smaller using density on 1 — 2 Mpc 
scales. Obiuc is also dependent on density in the sense that 
the blue peak appears broader for galaxies in denser envi- 
ronments, possibly more so at fainter luminosities. 

To demonstrate the reality of these trends, we show the 
actual colour distributions with overplotted fits in Figure [6] 
The strong luminosity dependence of /ibiuc, Cblue, /irod and 
(T re d and the density dependence of f rct j are immediately 
obvious. To better illustrate the differences, the values of 
/^biuc and p re d are marked as dashed vertical lines. Simi- 
larly the density dependence of <Tbi uc and er red can be seen 
by comparing the actual fit with the fit as it would be if 
a was fixed to the values computed for the lowest density 
bin (overplotted dotted lines). It is apparent, especially in 
the faintest, highest density bin (lower right) and for u — r 
near /ibiue and fi Te d th at a broader b l ue pea k and narrower 
red peak are required. iMasiedi et all l|2006h find that very 
nearby neighbours (within ~ 20") significantly influence the 
measured luminosity of galaxies: however model colours are 
likely to be much more robust. In any case, whilst Obiuc 
increases at high density, this is unlikely to be related to 
increased photometric errors: a le d should be as influenced 
but decreases, whilst the ~ 20% drop in the amplitude of 
the blue distribution (bottom right panel of Figure [6l com- 
pare solid and dashed lines) is too strong to be caused by 
such a rare population of very close neighbours. It must be 
noted that for such well populated samples there are small 
but mildly significant deviations from a perfect double gaus- 
sian distribution, visible as assymetries in the distribution 
or an excess of galaxies at intermediate colours. Whilst these 
deviations are potentially interesting, they are so small that 
they can and should be ignored in the context of this paper. 

The similar dependence of f ro( j on density across this lu- 
minosity range (i.e., no luminosit y dependence for A(f re d)) 
is in qualitative agreement with iBalogh et all |2004bl ). al- 
though they appear to measure a larger dynamic range for 
fred as a function of the surface density within the 5 th near- 
est neighbour, £5, at fixed luminosity. Binning more finely 
in density, we find that the increase in A(f re d) for galaxies 
in the highest density third is mainly contributed by a small 
fraction of galaxies (~ 15%) with £0,1 > 4.0Mpc~ 2 . It is 
also noticeable that the trends of A(/i rcd ) and A(^biue) with 
density are quantitatively comparable to the Balojjh et al.l 



|2004br i results, with more extreme values when more ex- 
treme (rare) environments are isolated. 

A detailed interpretation of our results will be left for 
Section [6] For now we note that Figure [3] suggests that 
suppression of star formation leads to the formation of red 
galaxies preferentially in environments characterized by a 
high local (< 1 Mpc scale) density. The importance of differ- 
ent scales will be addressed in Section f5. 21 The dependence 
of /i re d and <r re d on environment could relate to differences 
in either/both age or/and metallicity of the stellar pop- 
ulation, although intragalactic dust in galactic and group 
scale haloes can account for part of the eff ect (< 0.01 mag, 
iMenard et al.ll2009l : lM~cGee fc Baloghll2009l '). A plausible sce- 
nario is that less luminous galaxies in dense environments 
will typically be satellites, which may have experienced early 
suppression of their star formation. Such galaxies will have 
an older and purer underlying population than similar galax- 
ies in lower density environments, leading to redder colours 
and smaller scatter (i.e. smaller <r ro d). Meanwhile the redder, 
broader blue peak in dense environments suggests a more 
disturbed and (on average) less active recent star forma- 
tion history compared to those in less dense environments, 
although internal dust extinction may also play a role. 

5.2 Scale Dependence 

Figures [3] S] and [S] illustrate that whilst the colour distri- 
bution of galaxies depends upon density on < 1 Mpc scales, 
£0,1, it also depends (to a smaller extent) on a totally in- 
dependent measurement of density on larger scales, £1,2- 
However, although £1,2 is measured independently of £0,1, 
the two quantities are nonetheless strongly correlated. This 
is illustrated in Figure [7] in which we show the distribu- 
tion of —21.5 < M r < —20.0 galaxies on density computed 
on the different scales that we will consider in our analysis. 
Binning is logarithmic in density, with less than 50% density 
increase within each bin on either scale except the highest 
and lowest density bins, which have infinity and as upper 
and lower boundaries, respectively. The left panel examines 
the distribution in the £0,0.5 versus £0.5,1 plane, whilst suc- 
cessively larger scales are examined in the £0,1 versus £1,2 
and £0,2 versus £2,3 planes in the central and right panels 
respectively. 

The strong correlation between small and large scale 
densities at < 3 Mpc probes the non-linear regime of col- 
lapse which amplifies overdensities via infall of galaxies onto 
structures. For massive systems this correlation is further 
enhanced, particularly on small scales, by: 

• The presence of galaxies within the same virialized 
structure; 

• Projection effects: physical separation ^ projected sep- 
aration; 

• Redshift distortions (the "finger of God effect"): 
the true physical separation from neighbours within 
±1000 kms -1 will be a strong function of environment, with 
different contributions of galaxies inside and outside the lo- 
cal halo (in group catalogue construction, galaxies outside 
the local group halo are known as "interlopers"). However, 
a cylindrical depth of ±1000 km s -1 optimally red uces these 
effects with a fixed aperture (|Cooper et al.ll2005h . 

To study how the colour distribution of galaxies is af- 
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Figure 6. u — r colour distribution of galaxies (points and black solid line) as a function of luminosity and density, with double gaussian 
fits (red+blue peak=green line). Top to bottom: Decreasing luminosity M r = —21.5 to —20.0 in bins of 0.5mag width. Left to right: 
Undcrdcnsc to overdense, split evenly into bins by £o,i- Overplotted is the mean colour for each peak (vertical lines), indicating that 
both peaks move redwards going to higher density. Also overplotted (dotted lines) is the sum of two gaussians with equivalent parameters, 
except the widths ( (Tblue> "red) ar e fixed to the value computed for the lowest density bins (left panels). The deviation in the highest 
density bins illustrates the dependence of <r on density (strongest in the faintest bin, lower right). 



fected by density on different scales we perform the follow- 
ing analysis. First we bin the galaxies in the volume lim- 
ited sample according to their luminosity, in bins of 0.5 mag 
over the range 21.5 < M r < —20.0. Within each luminosity 
bin, galaxies are binned further according to the density of 
their environment computed in two spatial scales. Finally, 
the colour distribution of galaxies within each bin of lu- 
minosity and density on two scales is parameterized using 
the double gaussian model. Fitting of the characteristic pa- 
rameters is performed independently for each bin (provided 
that at least 25 galaxies are included) and residuals are 



computed to remove the luminosity dependence (A(f re d), 
A(ju r ed), A(o- r ed), A(^biue), A(cTbiue) )■ In order to improve 
the signal-to-noise of the parameter distributions as a func- 
tion of densities, we average the values of A in the three 
luminosity bins, weighted by the number of contributing 
galaxies. 

We illustrate that the fits thus produced are sufficiently 
good in Figure [8] To do so we must select a single luminosity 
bin, -20.5 < M r < -20.0. The top-left panel shows that the 
fraction of red galaxies, f ro d can clearly be traced as a func- 
tion of density computed on 0-0.5 Mpc vs 0.5-1 Mpc scales. 
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Figure 7. Number of 21.5 < M r < —20.0 galaxies binned in two independently measured scales of density. Left: on scales So, 0.5 an d 
So.5,i; Centre: on scales So,i an d Si ai Right: on scales £0,2 an< ^ ^2,3. Small and large scale densities are strongly correlated with 
each other. 
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Figure 8. Example double gaussian fits to the colour distribution of —20.5 ^ M r ^ —20.0 galaxies binned by density on scales So, 0.5 
and So. 5,1. Top left: The density dependence of the red galaxy fraction, f rcd for galaxies on scales So, 0.5 an£ i So. 5,1, as the left-hand 
panel of Figurc[9] except limited to galaxies in the —20.5 ^ M r ^ —20.0 luminosity bin. Other panels: Fits to the u — r colour distribution 
for galaxies in the bins indicated using black and white squares in the top-left panel. Panels refer to the bin in the equivalent position 
in the density-density plane. The colour distribution (points and black line) is in each case fit by a double-gaussian (green line) which 
is the combination of blue and red peaks (represented by blue and red lines). v/(iV + 2) errors are assigned to the counts in each bin of 
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To illustrate fitted populations at a variety of multiscale den- 
sities, and in particular including bins containing close to our 
minimum number of galaxies (25), we select the seven bins 
from this density-density plane as indicated by black and 
white squares. For each of these bins the full double gaus- 
sian fit to the it — r colour distribution is shown in the other 
panels. Most importantly, both peaks can be seen in each 
panel, and the fits are as good as the data allows. For bins 
containing few galaxies, especially those with few blue (or 
red) galaxies, the fit parameters are uncertain. Whilst this 
is formally accounted for in fitting errors (see Section [4]), 
errors on f ro( j are computed statistically sing the ICehrelsl 
(| 19861 ) approximation, combining Poissonian and Binomial 
statistical errors. For double gaussian fits to a large num- 
ber of galaxies this approximates the fitting error, whilst for 
smaller numbers the statistical error is larger. However, the 
scatter between different, independent bins in the density- 
density plane provides the best estimate of the "true" error: 
this is consistent with the statistical error, and is reduced 
onc e the three lumino sit y bin s are coadded. 

iKauffmann et all ( 2004 ) iBlanton et alj (120061 ) and 
iBlanton fc Berlindl ( 2007 ) have all claimed that galaxy prop- 
erties exhibit little dependence on density beyond f Mpc. 
Beyond this separation, the correlation function of blue or 
red galaxies is consistent within 5% of that expected if the 
type of galax y only responds to its ha lo mass and halocen- 
tric position (|Blanton fc BerlindllioOTl ). Our choice of pairs 
of small and large scales, whose effects are to be compared 
(namely 0-0.5 Mpc vs 0.5-f , 0-1 vs 1-2, 0-2 vs 2-3), is de- 
signed to provide complimentary information. In each case 
we have selected an "interior" (circular) annulus and a larger 
scale annulus which meet but do not overlap. This provides 
independently measured quantities, but means that any sig- 
nal on scales just beyond that measured by the small scale 
should be picked up on the large scale. To facilitate discus- 
sion of these scales, we introduce the division radius, rdi v 
which is equal to 0.5, 1 and 2 Mpc respectively for our three 
pairs of small and large scales. The Eo,i versus Ei,2 plane 
examines the influence of > 1 Mpc scales. Smaller scales 
(< 0.5 Mpc) trace local substructure. This is examined in 
the So, 0.5 versus So. 5.1 plane. Conversely, the Eq,2 and 
£2,3 plane probes the largest scales measured. If galaxies 
only correlate with the properties of their embedding halo, 
then one might expect large scale density only to matter 
where those scales trace galaxies within the same halo (di- 
vision radius, r^iv < 2r2oo, a typical halo diameter). This 
occurs at M ha i ~2x 10 13 M Q (5 x 10 12 M Q ,8 x 10 13 M Q ) 
for r div = IMpc (0.5 Mpc,2 Mpc). 

Figures O [TU] and QT] examine how the relative frac- 
tion of red galaxies ( A(f rc( j)), the parameters describing the 
red peak (A(p re d), A(cr rc d)) and the blue peak ( A(p;biue), 
A((Jbiu C )) depend upon small and large scale density for our 
selected planes of scale. Colour-coded bins in this plane are 
totally independent of each other. Overplotted contours are 
computed for a smoothed map of the data, using a gaussian 
smoothing kernel with a = 1.5pixels 0. Contours indicate 
the direction of change for each parameter in density-density 



space, and thus the importance of the two different scales. 
Smoothing reduces small scale noise, but note that local or 
abrupt features can be smoothed away. 

We also quantify the significance of any dependence on 
large scale density, at fixed small scale density. For each pa- 
rameter and set of scales, we compute the Spearman Rank 
Correlation Coefficient (p) of the correlation with large scale 
density, and the probability, Pl, of attaining the value of p 
given a null hypothesis in which the parameter depends only 
upon small scale density. The method used to compute Pl is 
described in Appendix[B] A value of Pl = 0.5 means no sig- 
nificance, while Pl = is a high significance anti-correlation 
and Pl = 1 is a high significance positive correlation 0. Sub- 
sections of the small scale density range are also tested so 
that effects of large scale density can be detected which are 
local in density space. The results are recorded in Table [T] 

A common feature of all the panels in Fig. [9l [lOl and ITTI 
is that contours, where they exist, are more closely aligned 
with the large scale axis than with the small scale one, in- 
dicating that the parameters react more strongly to smaller 
scales (even < 0.5 Mpc) than to the larger ones. 

Nonetheless, there is a highly significant positive cor- 
relation of A(f rc d) with Eo.5,1, at fixed £0,0.5- This must 
reflect a transformation process for galaxies which relates to 
the environment and is imprinted on scales beyond 0.5 Mpc. 
However parameters defining the red peak (A(/x re d) and 
A(cr ro d)) show no significant trend with Eo.5,1 at fixed 
So, 0.5, other than A(p rc d) at low density (£0,0.5 = — 
2 Mpc -2 ). A significant residual positive correlation does ex- 
ist for the blue peak position A(pbiuc) with Eo.5,1, though 
this is weaker than the trend with A(f re d). 

On > 1.0 Mpc scales there is no significant correla- 
tion of A(f r ed) with Ei, 2 at fixed Eo,i other than a ~ 2a 
level anti-correlation at low density ( Eq,i = 2 — 10 Mpc - ). 
A stronger anti-correlation exists between the position of 
the red peak A(p rc d) and Ei,2 at intermediate density 
(Eo,i = 1 — 10 Mpc -2 ). The blue peak position A(/ibiue) 
does correlate positively with Ei,2, though this is only sig- 
nificant at low density (Eo,i = — 2 Mpc -2 ). The lack of 
any positive correlation with densities above 1 Mpc scales, at 
least for A(f re d) and A(p rc d), and at higher densities typical 
of groups and clusters, indicates that there is no longer any 
direct impact of these environments on th ese galaxy proper- 
ties. This is consisten t with the results of Kauffmann et alj 
(|2004 ). lBlanton et al.l (|2006l 'l and lBlanton fc Berlindl (|2007t ). 
The anti-correlations are interesting however, and deserve 
further attention. 

At > 2 Mpc scales and at intermediate density ( Eo,2 = 
2— 10 Mpc -2 ), the anti-correlation of A(f re d) with E2.3 is 
even more significant. For roughly the same range of den- 
sity (E ,2 = 1 - 10 Mpc -2 ) we also see a highly significant 
anti-correlation of A(p rc d) with E2,3. In other words, red 
galaxies which have intermediate to large — 2 Mpc scale 
densities avoid regions overdense on 2 — 3 Mpc scales, and 
are otherwise bluer on average. Blue galaxies, on the other 
hand, exhibit a positive correlation of A(pbiuc) with E2,3, 
but only at low density (So, 2 = — 2 Mpc -2 ), as seen on 
1 — 2 Mpc scales. The width of the red peak A(<r r cd) is also 



6 The smoothing kernel is computed within a 5 X 5pixel box, with 
pixels beyond the data range assigned the same value as the edge 
pixel. 



Pl is computed with a numerical accuracy of 0.0003 
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Figure 9. The density dependence of the red galaxy fraction, A(f re( j), for 21.5 < M r < —20.0 galaxies binned in two scales of density. 
Left: on scales So, 0.5 an d So. 5,1; Centre: on scales Eo,i and Si 2; Right: on scales So, 2 and £2,3. Almost vertical contours illustrate 
that A(f rec j) depends mostly on small scale density. 
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Figure 10. The density dependence of the red peak parameters ^t rct j (above) and <r re( j (below), for 21.5 < M r < —20.0 galaxies binned 
in two scales of density. Left: on scales So, 0.5 and So. 5,1; Centre: on scales Sq,i and Ei,2; Right: on scales So, 2 and S2,3. 



anti-correlated (in this case the same direction as the trend 
with small scale density) with £2,3 at these densities. 



5.2.1 Summary of observational results: 

To summarize the way in which the galaxy colour distribu- 
tion traces density on different scales: 

• The fraction of red galaxies f re d increases dramatically 
with increasing local (< 0.5 Mpc scale) density. The proper- 
ties of both peaks also change: the red peak becomes red- 
der and narrower, and the blue peak becomes redder and 
broader. 

• The position and width of the red peak, /i rc d and er re d 
show no residual dependence on 0.5 — 1.0 Mpc scales. 

• The fraction of red galaxies and the position of the blue 
peak, f re d and /Ubiue, both show residual positive correlation 
with density on 0.5 — 1.0 Mpc scales. Other than some depen- 
dence of /ibiuc at very low densities, there are no significant 
positive correlations with densities computed on > 1.0 Mpc 
scales. 



• Nonetheless, at fixed smaller scales, there exist signifi- 
cant anti- correlations with large scale density (f re d and /i re d 
with Ei, 2, £2,3), typically at intermediate values of interior 
density. 

• Correlations with larger scale density exist at low values 
of small scale density (/i rc d with £0.5,1 at low £0,0.5, o" rc d 
with £2,3 at low £0,2, Mbiuo with £0.5,1, £1,2, £2,3 at low 
£0,0.5, £0,1, £0,2 respectively). This suggests that relatively 
isolated galaxies react in some ways to their environments 
on large scales. However this result might be explained by 
the effects of correlated noise on density measurements. We 
discuss this effect in Appendix [X] 



6 MULTISCALE DENSITY DEPENDENCES IN 
THE HALO MODEL FRAMEWORK 

6.1 Dependence of halo properties on density 

To aid our interpretation of these results, we cross-co r relate 
the sample with the groups catalogue of I Yang et al. I (|2007l . 
Y07) to examine how halo mass and group-centric radius 
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Table 1. For different choices of small and large scale density and a given range of small scale density (columns 1 and 2, units Mpc — 2 ), 
we provide the Spearman Rank Correlation CoefEecient p which measures the correlation between each parameter and the large scale 
density (columns 3 to 7). Pl, the probability that this coefficient is greater than would be expected in the absence of any intrinsic large 
scale correlation (i.e. any remaining correlation is purely a result of correlations between the parameter and small scale density, and 
between density on the two scales) is provided in columns 8 to 12. Computed using a monte-carlo technique, a value of Pl = 1 implies 
a highly significant positive correlation; Pl = implies a highly significant negative correlation and Pl = 0.5 implies no correlation. 
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depend upon our density measurements at different scales. 
We emphasize that our approach has been to compute an 
entirely model-independent, measurable quantity (density) 
and that this exercise thus serves as an examination of the 
halo model in an independent parameter space. 

The Y07 catalogue is constructed from the SDSS Data 
Release 4 (DR4) catalogue (via the New Y ork University 
Value Added Galaxy Catalogue, NYU-VAGC. lBlanton et all 



incomplete for our sample. A friends-of-friends al- 
gorithm is applied to the survey for group detection, and 
the total luminosity, L19.5, of galaxies is computed for each 
group down to an r-band magnitude of M r ,o.i C v = —19.5 
(note that in Y07 all luminosities are k- and evolution- 
corrected to 2 = 0.1). A halo mass is assigned to each group 
assuming a one to one relation in rank order of halo mass 
with 1/19.5 . This is in turn used to redefine group member- 
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ship, which is iterated until a stable solution is found. A halo 
is defined so long as a member galaxy with Mr-.o.iev < —19.5 
makes it into the sample. At the upper redshift boundary of 
our sample z = 0.08, this limit ensures that the group cata- 
log is complete to log 10 ( Mhaio/ M©) = 11.6 and includes all 
galaxies with M r < —20.434 (mag as defined in Section [5J 
corrected to z = rest- frame, no evolution correction). This 
sample, after removal of galaxies due to DR4 edge effects 
for which Mhaio could not be estimated, numbers 38733 
galaxies. Importantly, this catalogue provides halo masses 
for individual galaxies, as well as groups with more than 
one known member. 

Figures 1121 1131 and [14] illustrate how the halo masses 
Mhaio and halo-centric radii rh correlate with density com- 
puted in our three choices of small versus large scale density 
plane (rdi v = 0.5,1,2 Mpc). The left panels illustrate (from 
top to bottom) the dependence on density of the median halo 
mass, the fraction of galaxies in haloes of > 10 12 ' 5 Mq (both 
restricted to grid pixels with ^ 25 galaxies), and the me- 
dian halo-centric radius for galaxies in massive io 14 ~ 14 - 5 Mq 
haloes (for ^ 5 galaxies) . Each panel is also overplotted with 
contours, computed for a grid smoothed in exactly the same 
way as those in Figures 151 to 1111 

For direct visual comparison, the right hand panels re- 
produce the equivalent dependence of A(f re d), A(/x re d) and 
A(^biuo) on density, as already seen in Figures l9l to II II 

The median halo mass and the fraction of galaxies in 

> 10 12 ' 5 Mq haloes both depend mostly on the smallest scale 
densities. In more massive haloes the halo diameter exceeds 
r div = 0.5 Mpc (2r 2 oo ~ 0.4 Mpc for 10 12 5 M haloes) and 
thus So. 5,1 becomes more important, particularly at higher 
density where more massive haloes dominate. On Eo,i ver- 
sus Ei,2 scales (rdiv = IMpc), the contours are closer to 
vertical, and on Eo,2 versus E2,3 scales ( rdiv = 2 Mpc) both 
parameters anti-correlate with E2,3 at fixed Eo,2- This re- 
sults mainly from the way halo centric radius rh tracks den- 
sity on different scales. For io 14-14 - 5 Mq haloes the con- 
tours tracking the change in rh with density are almost 
vertical with rdiv = 0.5 Mpc, but moving to larger scales 
become diagonal for rdiv = 1 Mpc (larger median radius at 
smaller Eo,i but larger Ei^), and almost horizontal with 
rdiv = 2 Mpc (radius positively correlates with E2,3). The 
reason for this is fairly obvious: a galaxy at the centre of 
its halo will have high densities within annuli up to the size 
of the halo (r 20 o ~ 1.1,2.5 Mpc for M ha io = 10 14 ,10 14 5 Mq) 
and low densities on larger scales. A galaxy in the outskirts 
of this halo will experience the opposite, since the dense 
halo core is located within the larger scale annulus. Mov- 
ing to smaller halo masses than io 14-14 - 5 Mq, this pattern 
simply shifts to smaller scales. 

At the point of accretion, a galaxy moves from the in- 
falling halo (low Mhaio), to the more massive halo (high 
Mhaio and high rh), but the change in density will not be in- 
stantaneous. Therefore at high density on large scales there 
are not only more galaxies with large rh within the massive 
halo, but also more galaxies assigned instead to an infalling 
halo of lower mass. This can account for the anti-correlations 
of both median halo mass and the fraction of galaxies in 

> 10 12 ' 5 Mq haloes with large scale density. 

This effect is further amplified by the depth of the 
cylinder, required to account for redshift-space distortions 
(the "finger of God effect"). By integrating to ±1000 km s -1 



along the line of sight, the contribution of galaxies from out- 
side a massive halo is increased, including galaxies from the 
fore- and background. 

6.2 Comparison with the dependence of galaxy 
properties on density 

Now we make a qualitative comparison of the density depen- 
dence of galaxy properties with the density dependence of 
the assumed embedding halo to examine whether and how 
the halo model can explain the way the colour distribution 
of galaxies traces density on different scales. 

Both Mhaio and rh trace mostly < 0.5 Mpc scales in 
the Eo,o.5 versus Eo.s.i plane. Therefore any such depen- 
dence of galaxy properties can be explained by invoking ei- 
ther trends with halo mass or halo-centric radius. On larger 
scales (larger rdiv) this is no longer the case. The primary 
dependence on small scales (contours close to vertical) of 
A(f rc d), A(p re d) and A(p D iue) more closely resemble trends 
with halo mass than with radius. 

Indeed A(f re d) trends in particular are very well 
matched qualitatively to the trends seen for the fraction of 
galaxies in > 10 12 ' 5 M haloes, including the roughly equal 
spacing of contours. Trends on all scales are very compara- 
ble, but the anti-correlation with E2,3 at fixed Eo,2 is par- 
ticularly striking. The anti-correlation for haloes traces the 
accretion of galaxies onto the halo in the infall regions, and 
the reproduction of this trend for the fraction of red galaxies 
is intriguing. It provides strong circumstantial evidence that 
galaxies remain blue until they have experienced some trun- 
cation mechanism which is intimitely linked to the accre- 
tion time onto a massive halo. This is fully consistent with a 
simple model in which galaxies become red at some interval 
after their incorporation into a halo of mass > 10 12-13 Mq, 
which can also match t otal red fractions a s a function of 
halo mass and redshift (|McGee et al.l |2009| ) . However the 
qualitative dependence on density is not strongly sensitive 
to the choice of limiting halo mass: similar looking trends 
are produced with limits in the range ~ ]_q12.s-i3.s Mq, with 
contours becoming slightly more horizontal with larger lim- 
iting mass. 

[i rc d does not trace density on > 0.5 Mpc scales in 
the same way as the median halo mass or the fraction of 
galaxies in haloes above a given threshold. Instead the pos- 
itive correlations with density are all on local < 0.5 Mpc 
scales, and /i rc d anti-correlates with density on 1.0 — 2.0 Mpc 
scales. The lack of positive density dependence scales greater 
than 0.5 Mpc can be related to subhalo or radial effects, 
in that the colour may be imprinted at a time when the 
galaxy's < 0.5 Mpc scale subhalo was not part of the cur- 
rent > 0.5 Mpc scale main halo. Alternatively, if the colour 
of red galaxies is influenced by some combination of halo 
mass and halo-centric radius such that in massive haloes 
galaxies with larger rh are likely to be bluer than those in 
the centre, then it is also possible to qualitatively reproduce 
the observed trends. 

The colour of a red (assuming passive) galaxy can be 
driven by either age (luminosity weighted time since stars 
were formed), or metallicity of the stars. If age is the driv- 
ing factor then the subhalo hypothesis is hard to reconcile 
with observations: i.e. If galaxies all became red within such 
small subhaloes then we cannot explain the trend of f re d 
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with > 0.5 Mpc scale density. Halo-centric radius, on the 
other hand, should correlate with the time since a galaxy was 
accreted onto the halo, and provides a possible solution. If 
metallicity is the driving factor, the galaxy's environment at 
the time it was forming stars should be more important: i.e. 
the subhalo. Thus we suggest that where /i rc d is concerned, 
the relative importance of halo-centric radius versus local 
density is equivalent to the relative importance of age ver- 
sus metallicity. The anti-correlation of /i rc d with density on 
1 — 2 Mpc scales suggests that the average colour of red galax- 
ies is bluer outside more massive haloes (> 2 x 10 13 M ). 
This is consistent with both hypotheses (red galaxies be- 
longing to smaller infalling haloes are younger and metal 
poorer). 

We note that these trends describe only the influence 
of environment on red galaxies in the luminosity range 
—21.5 ^ M r ^ —20.0, and that the way in which both 
age and metallicity depend upo n their halo mass s eems 
to be a function of stellar mass (|Pasauali et alj|2010h . In- 
deed, the local density of red sequence galaxies correlates 
with b oth age and metallic ity, even at fixed luminosity and 
colour (|Cooper et alll2010h . whilst there is no clear correla- 
tion with larger scale (6Mpc) density (|Thomas et alJl2O10L 
who also show thee contribution of rejuvenated galaxies de- 
creases with density). Thus this complex picture has to be 
taken with the caveat that not all red galaxies are passive, 
and not all passive red galaxies have comparible star forma- 
tion histories. 

The blue peak colour /ibiue depends on — 1 Mpc 
scale densities in the same way as f re d, and therefore a 
pure halo mass dependence can be similarly invoked. /ibiuc 
probably traces the level of recent star formation in a 
typical star forming galaxy, and so this means that such 
galaxies experience partial truncation of their fuel sup- 
plies inside massive haloes, consistent with observations of 
truncate d gas disks or anem i c spirals in the Virgo clus- 
ter (e.g. Giovanelli fc Havnesl 19851; Koopmann fc Kennevl 



120041 ; iGavazzi et al.ll2006l ; IChung et al.ll2009l . and references 
therein). It is currently unclear whether this is necessarily a 
step on the path to total truncation and transition to a red 
colour. 

/ibiuc exhibits residual large scale dependence at low 
values of small scale density, no matter which radius rdi v is 
chosen to separate small from large scales. This dependence 
cannot be pure halo mass dependence, because the median 
halo mass has such residual dependence on large scale only if 
i'div = 0.5 Mpc. This might merely be a result of correlated 
noise in the density plane (see Appendex [A}> but there is 
no reason this should apply to certain parameters and not 
others. Otherwise the suggestion is that relatively isolated 
blue galaxies are redder when they live in relatively over- 
dense larger scale environments. On such scales (including 
> 2.0 Mpc) this might relate to an earlier formation time of 
a galaxy living in a large scale overdensity coupled with a de- 
clining star formation history, or to a more efficient cooling 
of hot gas onto galaxies in more isolated environments. 

We conclude that the halo model can provide an excel- 
lent qualitative description of the dependence of the param- 
eters governing the galaxy colour distribution, in particular 
frcd, Mrcd and /ibiuc, on simple annular density measure- 
ments on a range of scales from — 3 Mpc. The data is 
consistent with a simple picture in which the density depen- 



dence of the fraction of red galaxies f rc d and the peak colour 
of blue galaxies /^biuc is governed mainly by the fraction of 
galaxies assigned to haloes above some threshold halo mass 
(10 12,5 M© provides a reasonable but not unique description 
in this case). This suggests that the gas supply regulating 
star formation can be partially or completely suppressed as a 
function of halo mass. The halo model is particularly good 
at explaining the anti-correlation of f rc d with large scale 
density. This can be explained if galaxies not yet or recently 
incorporated into a relatively massive halo have not experi- 
enced the blue to red transition. The lack of positive corre- 
lations with density beyond 0.5 Mpc suggests that the red 
peak colour /i re d traces only local density within subhaloes 
(explicable if /i re d traces metallicity) or with additional de- 
pendence on halo-centric radius (since if suppression time is 
linked to the time a galaxy was accreted onto a halo, and 
this is correlated with radius, then a red galaxy will be older 
and thus also redder towards the halo centres). 



7 CONCLUSIONS 

We have developed a new multiscale method to character- 
ize the dependence of galaxy properties on environment. 
This provides a model-independent, rich parameter space 
in which to evaluate galaxy properties. 

In this paper we have examined the dependence of the 
u — r colour distribution for galaxies in the luminosity range 
—21.5 ^ M r ^ —20 on multiscale density. This is param- 
eterized using the double gaussian model, with parameters 
frcd (the fraction of red galaxies), /i re d and oved (the posi- 
tion and width of the red peak) and /ibiue and Obiue (the 
position and width of the blue peak). 

We confirm and extend known trends with small scale 
(< 0.5 - 1 Mpc) density: 

• Galaxies in denser environments are more likely to be 
red (frcd), implying an almost complete truncation of star 
formation. 

• If they are still blue (star forming), then they are also 
likely to be relatively redder in denser environments ( /ibiue) 
although with greater scatter ((Tblue), implying partial but 
varied truncation of gas supplies, possibly relating to trun- 
cated gas disks such as those observed in the Virgo cluster. 

• Red Galaxies are redder (/i rc d) and with less scatter 
( (Tred) in denser environments. This implies that they either 
form earlier and/or are more metal-rich, and that this for- 
mation path is also less varied than galaxies in less dense en- 
vironments. This depends upon the environment of a galaxy 
during its epoch of star formation. 

On larger scales: 

• No parameter correlates positively and significantly 
with density on scales > 1 Mpc, excluding trends seen for rel- 
atively isolated galaxies where correlated noise in the mea- 
surement of_d^nsity_jrnght be im p ortant. This confirm s the 
results oflKauffmarm et all |2004 ), lBl"anton et alj (^OOfj ) and 
iBlanton fc Berlindl (|2007l ). 



• Although all parameters correlate most strongly with 
density on < 0.5 Mpc scales, a residual positive correlation 
of frcd and /ibiuc with 0.5 — 1 Mpc scale density at fixed 
smaller scale density implies that total or partial truncation 
of star formation can relate to a galaxy's environment on 
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Figure 12. A comparison between the dependence of halo mass and halo-centric radius with density on scales So, 0.5 an d So. 5,1 (left) 
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panels on the left illustrate the behaviour of the median halo mass, the fraction of galaxies inside haloes > 10 12 ' 5 Mq, and the median 
halo-centric radius in haloes of lo 14 - 0-14 - 5 Mq, taken from the Y07 catalogue. 



these scales. A 1 Mpc diameter halo corresponds to a mass 
of ~ 2 x 10 13 M Q and a typical crossing time on the order of 
4Gyr, and one or both of these parameters must relate to a 
truncation mechanism which is active on these scales. 

• On > 2 Mpc scales f rc( j anti-correlates with density at 
fixed smaller scale density, and /i rc d anti-correlates with den- 
sity on > 1 Mpc scales. In other words galaxies are more 
likely to have ongoing star formation, or if not then to be 
relatively blue, when the large scale density is relatively high 



for the density on smaller scales (they live at projected dis- 
tances ~ 1 — 3 Mpc from a local overdensity) . 

We interpret these trends qualitatively using the halo 
model. To make this comparison we utilize the halo cata- 
logue of Y07 which is constructed using a friends of friends 
group finder, and under the assumptions that galaxies live in 
haloes, and that the mass of these haloes has a one to one 
relationship with the total galaxy light. The properties of 
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embedding haloes has been cross-correlated with our galaxy 
catalogue to examine how halo properties trace multiscale 
density. By comparing this dependence to that of galaxy 
colours, we find: 

• Qualitatively, the multiscale dependence of f rc( j and 
Mbiuc is remarkably comparable to that seen in the fraction 
of galaxies living in haloes above a given threshold mass in 
the range ~ 10 12 ' 5-13 8 M . This is consistent with a sce- 
nario in which galaxies can experience truncation of their 
star formation at some time after they are accreted onto a 
halo. 



• The halo model offers a simple and yet profound ex- 
planation for anti-correlations with larger (particularly > 
2 Mpc) scale density. For galaxies at large radial distances 
from the centre of massive haloes the halo core is included 
in measurements of density on large scales. Thus this region 
of parameter space traces the accretion region of haloes (ef- 
fectively marking the crossover from the one to the two halo 
term in the correlation function). Anti-correlations on these 
scales therefore trace the accretion region of haloes. That 
many such galaxies are still blue implies that the influence 
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of environment is felt at or after accretion onto a massive 
halo. 



is negated due to a combination of halo mass and radial 
trends. 



• Mrcd correlates positively with density only on < 
0.5 Mpc scales, implying that the red colour is set by physics 
within subhaloes of this size. In this case we suggest metal- 
licity effects must drive these trends, since we know that 
age must correlate with larger scale densities: i.e. galaxies 
become red due to the influence of 0.5 — 1 Mpc scales ( f re< j 
correlates positively with £0.5,1). However we cannot rule 
out age effects if the 0.5 — 1 Mpc scale dependence of /i re d 



Whilst we have restricted ourselves in this paper to 
qualitative comparisons within the context of the halo 
model, we advocate the use of multiscale density to quan- 
titatively examine physically motivated models (from sim- 
ple halo occupation models to semi-analytically populated 
haloes). This provides a uniquely rich and observationally 
motivated parameter space to examine the dependence of 
galaxy properties on environment. 
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APPENDIX A: POISSON SAMPLED 
NEIGHBOURS IN A CORRELATED DENSITY 
FIELD 

The number of neighbours brighter than a given luminosity 
limit, and within 1000 kms" 1 of a given galaxy, is clearly 
not going to provide the perfect measure of its environment, 
either in the sense of "that which best correlates with the 
galaxy's properties" or in the sense of "that which traces 
the underlying dark matter density field". In most senses 
this is unimportant: it merely provides a measurable tracer 
of environment at the position of each galaxy and is the best 
that we can do em p iricall y, without assuming a halo model. 

iBlanton et all (|2006l ) has suggested that a Poisson- 
sampled galaxy field can lead to correlated noise in the small 
versus large scale density plane and therein to false correla- 
tions with large scale density. Effectively, sampling of rare, 



quantized tracers such as bright galaxies leads to large er- 
rors sampled from a Poisson distribution, particularly on the 
smallest scales and the lowest densities where these errors 
are largest. Assuming some perfect, underlying measure of 
environment exists, measurements with large errors might 
have residuals which correlates with the large scale density 
measurement, which has a higher signal to noise. In reality 
this is highly complex, relating to our observing galaxies at 
random positions along their orbits; projection effects; the 
correlated measurements of density for neighbouring galax- 
ies and the less than 100% deterministic nature of galaxy 
formation which allows a galaxy to enter our luminosity se- 
lected sample. Modelling all these effects and their influence 
on the error vector in our density planes is well beyond the 
scope of this paper. 

Nonetheless, we do not believe that correlated noise sig- 
nificantly influences the departure from vertical contours, 
since we do have some plots with contours which appear al- 
most vertical at all densities, as confirmed using the statisti- 
cal test outlined in Appendix [B] If correlated noise has any 
influence, it is most likely at small density, where N/^/(N) 
is smaller, and thus we shall treat any possible influence of 
large scales at small density with caution. 



APPENDIX B: SIGNIFICANCE OF 
PARAMETER DEPENDENCE ON LARGE 
SCALES 

To test the significance of residual large scale dependencies, 
we setup a null hypothesis: that only smaller scale density 
matters. Under this assumption, the predicted influence of 
large scale density can be tested using bootstrap catalogues, 
resampled from the real one as follows: The galaxies in each 
bin of luminosity are ranked in order of their smaller scale 
density. For every galaxy in the sample, one of the twenty 
nearest in rank is selected to replace it, regardless of its large 
scale density. This effectively shuffles the large scale infor- 
mation whilst preserving the small scales (imposing vertical 
contours in the small versus large scale density plane). The 
usual fitting analysis is applied to the resulting bootstrap 
catalogue, and the residual parameter set ( A(f re d), A(fj, Te ^), 
A(a re d), A(^biuc), A(crbiuc)) is computed at each point in 
the 2D density grid. Three thousand such catalogues are 
created, resulting in three thousand such parameter sets for 
comparison with the true parameter set. 

The Spear man Rank Correlation Coefficient, p 
(|Spearmanlll904l ) tests the strength of correlation between 
the rank of one variable and another, and therefore does not 
assume a linear correlation between the parameters them- 
selves. Computing this coefficient for the correlation between 
each derived parameter and the large scale density for all 
bootstrap catalogues defines the expected distribution of p 
in the case that the null hypothesis were true (i.e. no large 
scale dependence). Therefore the significance that additional 
dependence on large scales exists (on top of that induced 
by correlations with small scale density, averaged over the 
whole 2D density grid) can be quantified as the frequency of 
coefficients determined from the set of bootstrap catalogues 
which are more extreme than that computed for the real 
data. We define Pl as the fraction of bootstrap samples for 
which the real data coefficient is greater than the bootstrap 
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one. A value of Pl close to one infers additional positive 
correlation on large scales, whilst values close to zero infer 
additional negative correlation. 



